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Abstract 

We study the transverse momentum distribution for a pair of heavy quarks pro- 
duced in hadron-hadron interactions. Predictions for the large transverse momentum 
region are based on exact order QCD perturbation theory. For the small transverse 
momentum region, we use techniques for all orders resummation of leading logarith- 
mic contributions associated with initial state soft gluon radiation. The combination 
provides the transverse momentum distribution of heavy quark pairs for all transverse 
momenta. Explicit results are presented for bb pair production at the Fermilab Teva- 
tron collider and for cc pair production at fixed target energies. 
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1 Introduction 

The distribution in the transverse momentum of a heavy quark pair, QQ, produced in 
hadron-hadron interactions is of interest for elucidating the underlying quantum chromo- 
dynamics (QCD), and its understanding is important in studies of B—B mixing and CP 
violation at hadronic facilities. Unlike the case of pairs of heavy quarks produced in e^e" 
annihilation, QQ pairs created in hadron hadron collisions are often not in a back-to-back 
configuration (even in a plane transverse to the beam direction). The net transverse mo- 
mentum of the pair measures the imbalance between the transverse momenta of the Q and 
the Q. In this paper, we examine the quantitative description in QCD perturbation theory 
of the expected imbalance. Predictions for the region of large ^j, are based on exact order 
QCD perturbation theory. For the regions of small and modest g^, we employ an all orders 
resummation of leading logarithmic contributions associated with the emission of soft gluons 
from the intial-state partons that participate in the hard scattering process. This calculation 
addresses a practical question for heavy quark tagging at hadron facilities: if a Q is tagged 
with a given transverse momentum, what distribution in transverse momentum should one 
expect to observe for the associated Q ? 

We consider the process hadron + hadron Q -\- Q -\- X . In the simplest parton 
model description, the underlying hard scattering process is parton -\- parton Q + Q. 
At this level of approximation, no bremsstrahlung gluons are radiated from initial-state or 
final-state partons. If one neglects intrinsic transverse momentum of initial state partons, 
q± is zero. Therefore, in the simplest parton model description, Q and Q pairs would be 
produced back-to-back in the transverse plane. However, gluon bremsstrahlung is important 
in QCD and in general generates non-zero qj_ . 

The single particle inclusive differential cross section for heavy quark production has 
been studied in detail at next-to-leading order in QCD. We cite in particular the calculation 
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of the complete first order corrections to the dominant QCD production channels 
and ^, and comparisons with data on inclusive b-quark production from the UA10 
and CDF collaborations0. We mention also the recent work on resummation of leading 
logarithmic contributions for the one particle inclusive cross section^]. In the single particle 
inclusive approach, the kinematical variables of the heavy quark's (or antiquark's) partner 
and of the final state light partons are integrated over with the attendant limitation that 
it is not possible to examine quark-antiquark correlations or the cross section differential 
in the transverse momentum qt of the QQ pair. Correlations have been studied at leading 



order 10 and at next-to- leading order] 10, 11 1. One may expect that next-to- leading order 
QCD should provide reliable expectations for the distribution in qt at large qt- At small 
qti the relatively large mass mq of the heavy quark Q justifies perturbation theory, but the 
presence of the two disparate scales, mq and g^, requires care. 

Perturbative QCD has been used for a successful description of transverse momentum 



distributions in massive lepton-pair production, the Drell-Yan process[0. There one studies 
the reaction hadron + hadron ^ e"*" + e~ + X, where the electron-positron pair is detected 
and its transverse momentum q± is measuredfl^, |l^ |l^, [16|. Important for the quantitative 



description of the transverse momentum distribution at modest values of is the resumma- 
tion of logarithmic contributions associated with emission of soft gluons in the initial state of 
the hard scattering process quark + antiquark — > e^ + e^ + X. We follow closely the analogy 
with the Drell-Yan case for the reaction hadron + hadron — ^ Q + Q + X and concentrate on 
the transverse momentum distribution of the QQ quark pair. Heavy quark pair production 
is, however, more involved than massive lepton-pair production. New complications arise 
from soft gluon emission from the final-state heavy quarks, effects that are absent in the 
Drell-Yan reaction. We will argue that the resummation technique for dealing with initial 
gluon radiation should still be applicable in our case. Soft gluon emission from final state 
heavy quarks has been studied in [|l^ . 
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For the process hadron + hadron ^ Q + Q + X , we choose to study the cross section 
differential in the variables and q"^; M is the invariant mass of the QQ pair. Knowing M 
and q± , we can judge how far the QQ system is away from the back-to-back configuration 
in the transverse plane. To be more specific, when M is near the mass threshold 2m of the 
QQ pair, the momenta of the Q and Q are close to zero in the center of mass frame of the 
QQ system. Only a small amount of q± then suffices to put the Q and Q in a configuration 
that is non-back-to-back in the transverse plane in the laboratory system of reference. On 
the other hand, at large M, the Q and Q have large relative momentum in the center of 
mass frame of the pair. At large M, large q± would be needed to produce a QQ pair that is 
not-back-to-back in the transverse plane in the laboratory system. 

At large and moderate values of q'j_ ~ (9(M^), the QQ pair production cross section 
can be computed perturbatively as 

al{ai + a2as + a-iol -\ ). (1) 



dM^dqi 

At any fixed order of as and qj_ ^ 0, the cross section is well behaved after the hard scattering 
cross section has been properly defined. At low q'j_ ^ 0, however, the convergence of the 
perturbative series deteriorates. For small q'j_ ^ 0, the dominant contributions (i.e. the 
leading logarithmic contributions) to Eq. (|l]) have the form 

^ '6i«,ln(^) + 62«2ij^3. ^2) 



dM^dql ql V ' ql ' ' ' q\ 

The convergence of the theory is therefore governed by asln^(M^/g^) instead of otg- The 
logarithms arise through emission of soft and collinear gluons. At sufficiently low g^, 
ln^(M^/g^) is large even when dg is small and any fixed order calculation breaks down. In 
order to obtain a reliable prediction, one must resum the leading contributions to all orders 
in ag. 
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The remainder of this paper is organized as follows. In Sec. 2 we present the pertur- 
bative calculation of the distribution using exact order QCD matrix elements. We 
describe how to obtain the asymptotic expression at — > from the exact order matrix 
elements. In Sec. 3 we discuss the formalism for resumming the initial soft and coUinear 
gluon contributions. We match the resummed result in the low-g^. region to the exact 0{al) 
result in the high-g^ region. Results and examples of the g^. distribution for specific hadronic 
reactions are given in Sec. 4. 

2 Perturbative Calculation 

We begin our discussion with the hadronic reaction in which a QQ pair is produced: 

p{K,) +p{K2) ^ Q{pi) + Q{p2) + X (3) 

In this expression, p and p denote the proton and antiproton respectively. The quantity X 
stands for all the final particle states that we sum over so that the above process is semi- 
inclusive with respect to the outgoing final particles. We use capital letters for the momenta 
of the proton and antiproton to distinguish them from those of the quarks, antiquaries and 
gluons. The corresponding partonic subprocess can be written as 

P^{h) +Pj{k2) ^ Q{P1)+Q{P2)+X, (4) 

where pij are the initial partons from the proton and the antiproton. The relationship 
between the momenta of the hadrons and partons is 

^1 = ^iKu k2 - 6i^2 (5) 

The four-momentum of the system made up of the QQ pair is 

q^^p'l+P^ (6) 
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We also define several useful invariants for our calculation: 

S = {Kr + K^f 
s = {h + hY = 66^ 
= = (P1+P2)' 

2 ^ 2q-Kiq-K2 _ ^2 ^ 2g ■ fci g ■ fca _ ^2 

We remark that is the collision energy of the proton and antiproton; M is the invariant 

mass of the heavy quark QQ pair; is the square of the transverse momentum of the pair, 

equal to the square of the vector sum of the individual transverse momenta of the Q and Q; 

Y is the rapidity of the QQ pair in the center of mass frame of the proton and antiproton; 

and y is the rapidity of the QQ pair in the center of mass frame of the interacting initial 
partons. 

The differential cross section for reaction (^), is expressed as a convolution of partonic 
cross sections and parton distribution functions: 



Here Tmin = iy + Q± + <l±y/S, and (P&ij^QQx/dM'^dq'^dy is the fixed-order reduced 

partonic cross section obtained by first calculating Feynman diagrams up to a given fixed- 
order in QCD perturbation theory and then implementing a renormalization scheme to 
remove any ultraviolet divergences. The soft divergences at — > are cancelled between 
the virtual diagrams and the bremsstrahlung diagrams. The collinear divergences (mass 
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singularities) at — * are absorbed into the definition of tlie parton distribution functions. 
Therefore the fixed-order reduced partonic cross section is well behaved for any value of 
g^, in particular qj_ = 0. (Note that the apparent divergence at qj_ = in Eq. (2) is 
cancelled by a 5-function that arises from virtual diagrams.) An advantage of calculating 
the order differential cross section in terms of is that q\ can be used as a cutoff variable 
for both the infrared and coUinear divergences in the bremsstrahlung diagram calculations. 
This means that at any finite value of q'j_, the order differential cross section can be 
calculated by evaluating the bremsstrahlung diagrams without explicit concern about soft 
gluon cancellation and factorization of the coUinear singularity. 

We now consider QQ pair production in QCD perturbation theory. It proceeds by the 
following two reactions in the Born approximation (order a^): 



Symbols inside the parentheses denote the momentum assignments for the partons. The 
Feynman graphs that contribute to the Born amplitude are shown in Fig. 1. 

The magnitude squared of the Born amplitude, averaged over initial colors and spins 




(9) 



and summed over final colors and spins, can be written as[^ 




qq^QQ 




F^QED 




(10) 



with 



A, 



■QED 



s s 



QED 
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The total partonic cross sections are|T^ 



with 



:ii) 



ti and ui are defined as 
It is easy to verify that 



(12) 



Qq,99^QQ 
dM^dqldy 



(13) 



4^(M^) 



M2 Ar2 _ 1 



:i + 4 



M2 



Inx 



+ 



M2 Ar2_i 



, 1 5 . 4m^ , 
I -ID In X 

^ 3 3M2^^ 



(14) 



(15) 



/3 = a/1 -4mVM2, X 



1 + /3' 



(16) 



The Ca and Ci;' are the Casimir invariants for the adjoint and the fundamental representation 
of SU{N). For the particular case of SU{3), 



iV = 3, Ca = N, Cf 



2N 



(17) 



Observe that the differential cross section Eq. (|l^) is proportional to S{q'j_), because in the 
Born approximation the heavy quark pairs are produced with = 0. The total partonic 



cross sections ^^^''(M^) and a^^J{M'^) are inversely proportional to M^. This dependence can 
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be understood because the lowest order graphs either have only s-channel poles or fermion 
exchange lines. Most of the heavy quark pairs are therefore expected to be produced near 
threshold, = 4m^. 

At order the contributing partonic subprocesses include gluon bremsstrahlung dia- 
grams and (anti)quark-gluon scattering diagrams: 

q{h)+q{k2) ^ 9{h) + Q{pi) + Q{p2) 
g{ki)+g{k2) ^ g{h) + Q{j>i) + Qij>2) 
9{h) + q{q){k2) ^ q{q){k:,) + Q{pi) + Q{p2) 



plus virtual correction diagrams ||2CI|| to the lowest order processes Eq.(P). In Fig. 2 we show 
some examples of the gluon bremsstrahlung diagrams that contribute at order a^. 

In the order processes (p^Sj), q\ is no longer constrained to be zero; a spectrum of 
values of q\ will be produced. The q]_^ Q limit means that the parton with momentum k^ in 
the processes (0), is either soft and/or coUinear to one of the initial or final state partons. We 
note that in the qq channel there are diagrams with a gluon emitted from an initial quark or 
antiquark line; these look diagrammatically exactly the same as corresponding graphs in the 
Drell-Yan reaction. In addition there are diagrams with a gluon emitted from a final heavy 
quark or antiquark line; these are absent in the Drell-Yan reaction. The initial gluon emission 
diagrams can be soft and/or collinear divergent at — > 0. The final gluon emission diagrams 
can have only a soft divergence at — because the final state quarks are massive. The 
same statements can be made for the gluon gluon initiated subprocesses which are usually 
the dominant processes for QQ pair production. The quark-gluon scattering diagrams have 
only a collinear divergence from light quark emission. The square of the matrix element for 
the order processes Eq. (|T8|) has been published]^. We will make use of those results in 
our calculation. 



To calculate the partonic differential cross section, daf^_^QQ^/ dM'^dq\dy, at order 
we must integrate over variables which are independent of M^, q\ and y for the gluon 
bremsstrahlung and quark-gluon scattering processes. We choose the QQ center of mass 
frame in which 

= iPi+P2r = (M, 0, 0, 0) 



Pi = (M/2, cjo sin 6'i sin 6'2, cuq sin6'i cos6'2, co'ocos6' 



P2 = (M/2, — tuo sin sin 6*2, — tuo sin 6'i cos6'2, — u;o cos^^ 



(19) 



with 



uo = ^v'l-4mVM2. (20) 
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We obtain 



<^QO..(p,rt) = ] .., ^?^f(l-^)(l-?) ''Z"^' 



dM'^dqldy'-' ' 128 TT^ sM2 ^1 ^2 V ^2 l + gi/M2^ 



X 



The integral / £'p\ can be expressed explicitly as 



'^''''^((^-Pi)^--^) WUo.- (21) 



/^K(^-^i)'-^') I^I'.^qq;^ = - AmyM^J^^de, sine, j^e, \m\ 



■2 



(22) 



and |Af I .^.^ggj^. is square of the matrix element for the order processes (|T8D, averaged over 
spin and color. 

It is possible to perform the angular integrals analytically. However, the squared matrix 
elements for the QQ production process are lengthy and the integrals are difficult, making 
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the calculation formidable. Fortunately, for any finite value of q\_ the order differential 
cross section Eq. ( pT]) is free from any soft and coUinear divergences so that the angular 
integrals Eq. ( p2D can be dealt with numerically. Then the differential cross section for 
the proton antiproton reaction Eq. can be calculated by convoluting the partonic cross 
section Eq. ( pT]) with the parton distribution functions according to Eq. (^). 



A quantity of interest and importance is the average of the square of the transverse 



momentum at a fixed invariant mass M. It can be calculated asf22 



< ql{M) > 



da, 



pp^QQX 



, 2 2 d'^'^pp^QQX 

"^^^ dM^dql 



E /i) /i) / dql ql .Jf.^f , 

i,j't^,q,g-^Mys JMyUrS) '''^ Jo dM^dqi 



where 



(23) 



d(T. 



pp^QQX 
dM^ 



da 



(0) 



(1) 



ij—fQQ _|_ ij^QQX 



dM^ 



dM^ 



ijt^^qjMyS JA'P/i^iS) dM^ 



(24) 



Symbol /C denotes the familiar K-factor, not necessarily a constant in this case. The 
transverse momentum averaged over all values of invariant mass can also be calculated in 
the same way. 

As remarked in the Introduction, the convergence of the perturbation series deteriorates 
in the region qj_ <C M^. For predictions of improved reliability in that domain, one must try 
to sum leading contributions from all orders in a^. We will use the procedure for resumming 
contributions from initial state soft and collinear gluons developed by Collins and Soper. To 
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begin, we must first extract the leading contributions in the region of small q\ <^ from 
the cross section. This is done by calculating the angular integrals in Eq. (p2D analytically. 
Fortunately, in the limit g^/M^ — ^ the squared matrix element simplifies substantially and 
the angular integrals can be calculated analytically without too much difficulty]^. 

We take the limit — of the exact expressions for the square of the order matrix 
elements and do the integral under the same limit. The calculation is accomplished in two 
steps. First, we calculate the soft gluon contributions by setting the gluon momentum ^ Q 
everywhere in the matrix elements except in the denominators that are singular as /cs ^ 0. 
The soft gluon matrix elements have also been derived in the literature (see Eqns.(5.1-5) in 
Ref.pi and Eqns. (2.24-26) in Ref.Q]). In the limit 0, the 2^3 kinematics can also 

be approximated by 2 — > 2 kinematics, which we can implement effectively by replacing the 
5— function in Eq. (^TD by 



The singularities in the denominator can be replaced by 

=^ il. (26) 

Second, the hard coUinear contributions can be calculated if we replace the first 5— function 
in Eq. dll) by 

V ^2 1 + (l-X2/6)+^(l-Xi/6)+ ^ ^ 

The "+" prescription in the above equation, defined as 

^ fix) , fix) - /(I) , , , 

■'^ ' dx= ^ ^rfx, (28) 



{l — x)+ Jo 1 — x 
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ensures that there is no double counting in the phase space region where the soft and colhnear 
divergences overlap. Both replacements, Eq. ( ^Sf ) and Eq. (pTl) , imply significant simplifi- 
cation for our calculation. Their origin has been well demonstrated, e.g. in Eqns.(2.11-13) 
of Ref. [^. After we perform the angular integrals and convolution with the parton dis- 
tribution functions we obtain the asymptotic expression for the differential cross section for 
process (0) in the form 



^'^^^"-«^^:(asym) 



dM^ dqldY 



S 271 ql 



+ 



fi/p{Xl) (Pj^h ® fb/p) {X2) + [Pi^a ® fa/p) {Xl)fj/p{x2) 



(29) 

The symbol ® denotes a convolution of the parton distribution function / and Altarelli-Parisi 
splitting function P, defined by 



f®p){x)^ f'mp{-)^. 

' ^y^ y 



(30) 



The functions aI]'' come purely from initial state soft and collinear gluon radiation: 



(31) 



The functions iS-J"* can be split into two parts, in terms of initial and final state gluon 
radiation: 



(32) 



with 



r{1) 



-3Cf, = 2Cp + 2C^^(l-2m7M 



Inx 



2 /i\/r2\ 



(3 
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/5 



r(i) _ n r(^) 

^q{q)g,I - ^' ^q(q)9,F 



(33) 



Quantities x and /5 have been defined in Eq. ([T6|). Our coefficient I3q of the /3-function is 
normahzed to 

/3o = ^ ^. (34) 

To verify our calculation we have checked that the 1/e pole terms generated by di- 

mensional regularization from the terms not proportional to the Altarelli-Parisi splitting 
functions in the expression Eq. (^9]) cancel the infrared pole terms from virtual diagram 
calculations IPI, ffl. 



The expressions for the initial state gluon radiation terms for the qq channel, A, 



(1) 
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and -Bg^^/, are exactly the same as those in the Drell-Yan reaction. For initial state gluon 
radiation in the gg channel, our expressions A^^^ and Bfj also agree with those found by 



Catani, D'Emilio, and Trentadue|p5|| and by Kauffmanp3| for production of a color singlet 
state from gluon-gluon fusion. The agreement in both cases follows the expectation that 
the initial state soft gluon radiation does not depend on the type of hard process under 



consideration [^, It indicates that the effects of initial state gluon radiation can be 



resummed to all orders of as for QQ pair production, as in the cases of the Drell-Yan 
reaction [13, 14] or Higgs production through gluon fusion]^, 25, 23 1. In the next section we 



will attempt the resummation of the initial soft and collinear gluon contributions using the 
formalism developed for the Drell-Yan process [^, p9[] . 

The expressions for ffnal state soft gluon radiation, B'^^^p and B^^Jp, are the same for 
the qq and gg channels. This can be understood since the ffnal state soft gluon radiation 
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occurs from the final state (anti)heavy quark lines in both the qq and gg channels. Soft 



gluon emission from final state heavy quarks has been studied in Ref. ||17||. 



3 Resummation 



The technique for resumming contributions from initial state soft and collinear gluons was de- 
veloped by Collins and Soper|2^, and it has been applied to massive lepton-pair pr o duct ion fl^, 



single vector boson production |3y, Higgs boson production!^, ^ |2^, and ZZ-pair 
production [^]. In our case, the appropriate expression analogous to that of Collins and 
Soper is 



(iM2 dql dY 



resum 



alf{M') ^ dbbJoibq^) W,,{M,b) , (35) 



where Jo{x) is the zeroth order Bessel function. The function Wij{M,b) sums all the loga- 
rithmic terms of the form a^ln™' {M'^b'^) with 1 < m < 2n in the impact parameter b space. 
The all orders structure of W is given by the functional form 



W,j{M,b) 



clM^ dq^ 



exp<^ ; 9 



ln(^^^) MaM')) + B,,{aM')) 



(J2 (J2 
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(36) 



The parameters Ci, C2 and C3 are somewhat arbitrary. They are associated with the choices 
of renormalization and factorization scales in a fixed order perturbative calculation. We use 
the standard choices 



C2 = l; 



(37) 
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where 7^; is Euler's constant. The symbol (S> denotes a convolution, defined in Eq. (30). In 
the limit q^/M 0, the parton momentum fractions are 



Y -Y 



The functions A, B, and C{x) may be expanded in a perturbation series in a^: 



= E4;^(^j ; (39) 

The reason for the bar over the j in the expression for Cij{x, as) is that the fiavors of i and 
j must be the same in the case of qq. 

We work to first order in the expansions of A and B which corresponds to summing 
the first two powers of ln(M^/g5.) at every order in ag, i.e, the double-leading logarithm 
approximation. The A^^^^s and -Bj-!*'*'s depend implicitly on the choices of Ci^2,3- The simplest 



forms result from the choices in Eq. (|37D. The coefficients a\^^ and B^j^ in Eq. ( |39D can 
be obtained by formally expanding Eq. (^) in a series in as and then comparing with 
the asymptotic perturbative calculation from our previous section. For initial state gluon 

radiation, the expressions for A,-]'* in Eq. (p9| ) are the A\f of Eq. (0), and the B^j^ are the 

B^^jj in Eq. (|33|). For C, we make the simplifying choice C^j^ = — x) since we are 

working in perturbation theory to the first non-trivial order in as for large q±. (Note that 
crf°^ in Eq. (35) is proportional to a^.) Our neglect of C^j^ will only affect the normalization 
at gj_ = to 0{as) and the distribution for g_|_ 7^ to 0{as)- These statements imply that 
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our calculation includes resummation of the double-leading logarithms to all orders in 
but the total integrated cross section is accurate only to second order in a^. 



We comment that the validity of the resummation formalism p9|] was demonstrated for 
the Drell-Yan reaction and for W, Z production where there is no final state gluon radiation. 
We are making the reasonable assumption here that the same formalism is valid for dealing 
with the effects of initial state soft and coUinear gluon radiation in the case of heavy quark 
pair production. Resummation of soft gluon emission from the final state heavy quarks has 



been studied in Ref. [17 



The gluon resummation formula, Eq. (pS]), provides the cross section in the region of 
small g^; for the high-g^ region we use the exact 0{al) perturbative calculation. We will 
join the results for the low q± and high q± regions using a matching procedure employed 
previously p3| . 



da da 

(match) = 2 (pert) 



dM^ dql dy ^ ' c/M2 dq\ dy 



The function 



da da 

(resum) - (asym) 



dM^dqldy' ' dM^dqldy 



.(40) 



f(q±/M) = — - (41) 

■'^^ ' ' l + (3gx/M)4 ^ ^ 



serves to switch smoothly from the matched formula to the perturbative formula p3[. For 
details of the matching procedure we refer to papers by Arnold and Kauffman and subsequent 
publications. 

Before presenting numerical evaluations, we end this section with a few remarks. As 
discussed by Parisi and Petronzio[Q, the resummed expression Eq. (|36D is ill-defined when 
h > Ij Kqcd because confinement sets in and blows up. Procedures have been proposed 



in the literature 1 28, 29, 14 to deal with this difficulty and parameterize non-perturbative 
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effects. In this paper we follow the method used by Collins and Soperp^ and by Davies and 



collaborators IjT^. We replace W{b) in Eq. (^) by 

W{b) ^W{b,)e-^-^^^^; (42) 

(43) 



1 + 62 /52 



Large values of b are thereby cut-off at some bmax', exp(— S'„p(6)) parameterizes the large-6 
dependence due to nonperturbative physics. In principle, exp(— S'„p(6)) can be measured, 
but in practice one can approximate the function with a simple Gaussian parametrization, 

Snp{b) = b^igi + 92 HbmaxM/2)]. (44) 

According to Davies and collaborators H 



gi = 0.15GeV\ g2 = OAGeV\ b^ax = i2GeVy\ (45) 

The values of gi and g2 are obtained by fitting massive lepton-pair production data at 

y/S = 27 and 62 GeV. There is no strong reason to believe that the contribution from non- 
perturbative intrinsic transverse momentum should be identical for subprocesses initiated 
by gluon-gluon scattering, as in our case, and quark-antiquark scattering, as in massive 
lepton-pair production. When substantial samples of data become available on cc and bb 
production, it should be possible to refine the choices made here. Particularly informative 
in this respect will be data on the azimuthal angle (0) dependence. The extent to which the 
quark and antiquark are produced with near vr is particularly sensitive to the net transverse 
momentum imparted to the quark-antiquark pair|lO, [Tl|. 



4 Results and Discussion 



In this section we present and discuss some phenomenological applications of our analysis. 
We use HMRS parton distribution functions [^] and the one- loop corrected formula for 
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the running coupling constant as(/i) with A4 = 190 MeV. The factorization scale and the 
renormalization scale are chosen to be the same as the invariant mass of the heavy quark 
QQ pair, i.e. = M, unless stated otherwise. 

In Fig. 3 we show the lowest order result for the distribution in invariant mass of a bb 

pair produced in a pp collision at the Fermilab collider energy = 1.8 TeV. The bottom 
quark mass is chosen as rrih = 4.75 GeV. The distribution peaks at a value of M a few GeV 
above the bb pair mass threshold of 2mh. It then decreases quickly as M increases. This 
implies that most b's and/or b's are produced at the Fermilab collider with small momentum 
in the bb center of mass frame. Next-to-leading order QCD contributions change the overall 
normalization of this curve, but more important for the sake of our present discussion, they 
should not change the shape of the lowest order curve except very near threshold or far above 
threshold [g. 

In Fig. 4, the average quantity < g^(M) > /C is plotted as a function of the pair 
invariant mass M; the factor /C was discussed in Sec. 2. The quantity < q±{M) > /C is 
proportional to the square of the average transverse momentum of the bb pair. Its value is 
about 80 GeV^ near threshold and rises linearly with M in the range shown in the figure. 
At the fixed rapidity value Y = 0, the function has the same shape and magnitude as in 
Fig. 4, understandable because the bb pairs are produced centrally. (An approximately linear 
rise with M of the average transverse momentum in the region r = A/P/S <^ 0.05 would 
be expected from simple dimensional arguments. The growth of the square in Fig. 4 is less 
rapid than quadratic.) 

An important inference may be drawn from Figs. 3 and 4. We may estimate from Fig. 3 
that the average invariant mass < M >= 15 GeV. Glancing at Fig. 4 we notice that near 
and above the pair mass threshold < q'j_{M) >/C is larger than 80 GeV^. Using /C ~ 2. 4 (see 
ref. m, H), we deduce that bb pairs produced at the Fermilab collider are expected to have 
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an average transverse momentum about 5 GeV. As will be discussed below, after integrating 
over all values of M, we find that the square of the average transverse momentum < > 
is 36.0 GeV'^ in the purely perturbative order approximation. Taking the square-root, 
we deduce < q± > rms'^ 6 GcV . This value is comparable to, and slightly larger than, the 
typical momentum of an individual 6 or 6 in the bb center of mass frame. Correspondingly, 
a significant fraction (about a quarter to half depending upon how one defines back-to- 
back configuration) of bb pairs at the Fermilab collider are expected to be produced in a 
configuration that is not back-to-back in the transverse plane. 

The differential cross section da/dMdq± is presented in Fig. 5 for three fixed values 
of mass, M = 15, 25 and 50 GeV. The dashed curves are from our fixed order per- 
turbative calculation. They are most applicable in the region q± ~ 0{M) where there is 
essentially only one hard scale in the problem. The fixed order results become inapplica- 
ble if qj_ << 0{M) where, as discussed earlier, the effects of soft gluon contributions must 
be incorporated in order to obtain a more reliable result. The dot-dashed lines show the 
asymptotic results, Eq. (29), obtained from the fixed-order results in the limit q± 0. At 
small q± the asymptotic results agree with the perturbative results, as expected. At larger 
q±, the asymptotic results manifest unphysical characteristics that can be traced to the fact 

that the functions A^f and B^j^ have opposite signs (c.f. Eqns. (31) and (33)). Accordingly, 
we take the asymptotic results at face value only at small q±. The cross sections obtained 
from resummation of the effects of initial soft gluon radiation are shown by the dotted lines. 

The resummed results are not expected to follow the asymptotic results because only 
the initial-state soft gluon radiation is included in our resummed formalism. The resummed 
and the purely perturbative order a'l curves nearly coincide for q± about 5 GeV and greater, 
as might be expected since the mass of the bottom quark is the relevant physical scale at small 
q±. Owing to the effects of resummation, the shapes of the two curves differ significantly for 
q± less than 5 GeV. 
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The solid lines in Fig. 5 present our final matched results, obtained from Eq. 
The matched results agree with the resummed results in the region of small q± and with the 
perturbative results at large q±. These three figures demonstrate how the final resummed 
and matched results differ from the perturbative results in the small g_L region. They also 
show that the simple matching procedure seems to work adequately in our case|^ . 

To preclude confusion, we stress that da/dMdq± is presented in Fig. 5; thus, the 
vanishing of the resummed curves as q± goes to zero has a kinematic origin. The divergence 
apparent in Eq. (2) is not present in the resummed calculation. 

In Fig. 6(a), the perturbative, asymptotic, and the resummed results from Fig. 5(a) 
are replotted as da/dM/dq'j_ versus qj_. This figure illustrates the behavior of the cross 
section at small q± in a different way, without the phase space factor of q± that is present 
in da/dM/dq± shown in Fig. 5(a). The same results are plotted again in Fig. 6(b) but 
with the switching function Eq. (41) included as a multiplicative factor in the asymptotic 
and resummed results. A comparison of Fig. 6(a) and 6(b) demonstrates the effects of the 
switching function included in the matching formula Eq. (40). At small q±, the switching 
function is close to 1 and does not modify the asymptotic and resummed results. At large 
q\ ~ 50 GeV^, the switching function suppresses the asymptotic and resummed results by 
almost a factor of 10. 

In Fig. 7 our final results are shown for the distribution in the square of the transverse 
momentum of the bb pair, da/dq\. Here we have integrated over M. As remarked at the 
start of this section, for consistency, all curves, including the purely perturbative order 
curve (dashed line), are computed with the one- loop evolved form for 

For g_L 7^ 0, the differences between the matched (solid curve) and fixed-order 
(dashed) results in Fig. 7 are formally of order and higher(except that our integrated cross 
section is valid only to 0{a1), as explained in Sec. 2). These differences will affect, among 
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other observables, the predicted average transverse momentum of the bb pair. Computations 
of the integrals over all q± of the product of qj_ times da/dq"^ are straightforward for both the 
solid and dashed curves in Fig. 7 since this product is finite and well behaved as q± approaches 
zero in both cases. The integral of da/dq]_ itself is straightforward for the matched case, 
where there is no divergence as q±_ approaches zero, but is more involved in the purely 
case. In the case, full account must be taken of virtual diagrams that contribute at q]_ = 
ID, H, If. Carrying out the computations, we obtain < q'j_> = 36.0 GeV^ in the purely 
perturbative order approximation and < q\ > = 66.7 GeV^ for our matched case. We 
have checked that use of the two-loop evolved form for as changes the purely perturbative 
order value of < > by less than 1 GeV^. 

Using the numbers in the paragraph above, we note that < > is increased by 
30.7 GeV^ as a result of soft-gluon resummation and matching. Taking square- roots, we 
find < q± >rms^ 8.2GeV in the matched case, to be compared to ~ 6 GeV in the purely 
perturbative order case. It may seem remarkable that an additional ~ 2 GeV is associated 
with soft-gluon resummation. It would be useful to be able to compare this predicted increase 
with that expected for < > in massive lepton pair production (the Drell-Yan process) at 
Tevatron energies, at massive lepton pair masses in the vicinity of 10 to 20 GeV, comparable 
to those for bb pair production. However, to our knowledge, no calculations have been 
published for the resummed q± distribution for the Drell-Yan process at such masses at 
Tevatron energies. Calculations at lower energies P5|, |T^ and/or higher masses PU| pT| may 



not be a useful guide since the distribution should broaden with energy at fixed mass. 
We remark, however, that the q± distribution should be significantly broader for bb pair 
production because the coefficients in the resummation expression, Eq. (36), are much larger: 

^« = 2Ga vs. 4^ = 2Gf. 

It would be interesting to compare the solid line in Fig. 7 with data. We comment 
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that our result is for production of a pair of b and b quarks, not a pair of B and B mesons. 
To compare with measurements of the transverse momentum of a pair of B and B mesons, 
one must include effects associated with fragmentation of the b and b quarks, and try to 
estimate effects associated with final state gluon radiation. In a Monte Carlo simulation 
of the single 6 or 6 inclusive spectrum, Kuebel and collaborators showed that the effects of 
final-state and initial state gluon radiation tended to compensate in some instances, with 
final-state radiation tending to soften the spectrum and initial-state radiation broadening the 



distribution |]IOl. To obtain the b quark inclusive cross section, the UAl collaboration fg] and 



the CDF collaboration |]36| use a Monte Carlo procedure to take into account the fragmenta- 



tion effects of b quarks into observed B mesons or single leptons before they compared their 
data with theoretical predictions. A similar procedure is required before comparison can be 
made of data with our theoretical prediction of the pair transverse momentum distribution. 

In Fig. 8, we show a distribution in the square of the transverse momentum for produc- 
tion of a pair of charm quarks, cc pair production in fixed-target proton beam experiment 
with Eheam = 800 GcV. In the case of charm, the average transverse momentum is so small 
that non-perturbative physics may dominate in the region of small q±. We will not dwell here 
on applications of the resummation method to charm pair production. In another paper, we 
plan to present further phenomenological applications for bb pair production. 

In this paper, we have focussed on the distribution in transverse momentum of a pair 
of heavy quarks. In carrying out our calculation, we integrated over (angular) variables in 
the QQ rest frame. Correspondingly, certain limitations must be accepted. In the context of 
our calculation, we are not able to describe the fully differential distribution in the momenta 
of the Q and Q separately, notably the azimuthal angle (0 ) dependence in the transverse 
plane nor correlations in rapidity. The limitation on the description of rapidity correlations 
appears insignificant since these differ little at leading and next-to-leading order [|TI| . On 



the other hand, the azimuthal angle dependence is sensitive to the net transverse momentum 
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imparted to the QQ pair [[TT], |T^, and it would be valuable to develop our approach further 



in order to examine the influence of soft-gluon resummation on the (f) distribution. 

A calculation of the fully exclusive parton cross section for QQ pair production at or- 
der has been published, along with examples of distributions at collider and fixed-target 
energies [|ll|]. As in the calculation reported here, those results include the lowest order a1 
cross section for production of a QQ pair, the order virtual corrections to the lowest order 
cross section, and the order cross section for production of a QQ pair along with a light 
parton. That calculation does not include the effects of soft gluon resummation presented 



in this paper. On the other hand, included in Ref. [pH]] is an exploration of certain other 
effects that go beyond the pure QCD calculation. The parton shower Monte Carlo pro- 
gram HERWIG was used to simulate the effects on the distribution of finite intrinsic 
transverse momentum of the initial partons. Substantial broadening of the distribution 
was observed, tantamount to that one would expect if the incident partons carried an intrin- 
sic transverse momentum of about 1.7 GeV. While such a large intrinsic contribution was 



questioned in Ref. [|TT| as perhaps an unreliable artifact of the parton shower algorithm, the 
notable influence of the added transverse momentum on the distribution underscores the 
importance of the type of study carried out in the present paper. We recall that our soft 
gluon resummation introduces substantial additional < g_L >. As a step in the direction of a 
full investigation of the influence of soft-gluon resummation on the distribution, it might 
be possible to incorporate the matched distributions shown in Figs. 7 and 8 into a modified 
order event generator. 

In summary, we have studied the distribution in the transverse momentum of a pair of 
heavy quarks produced in hadronic reactions. For large the order perturbative result 
should be applicable. In the region of small g^, we argued that resummation techniques 
developed in the study of the Drell-Yan reaction should apply for initial state soft gluon 
radiation. Use of the resummation method, plus a matching of results in the small and large 
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q± regions, permits an improved prediction of the full q± spectrum. Numerical results for 
the region of small q±_ region were presented for bb pair production at the Fermilab collider 
and for cc pair production in fixed target experiments. 
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Figure Captions 

Fig. 1. Feynman diagrams at order a^. 

Fig. 2. Examples of order Feynman diagrams with a gluon emitted from an initial quark 
or gluon line or from final heavy quark line. 

Fig. 3. bb pair invariant mass distribution computed from lowest-order QCD processes. 

Fig. 4. Average of the square of the transverse momentum of a bb pair, multiplied by the 
factor K, as a function of the pair invariant mass M. This curve is obtained from the 
purely perturbative order calculation. 

Fig. 5. bb pair transverse momentum distributions for three values of invariant mass, M = 
15, 25, and 50 GeV. We show our fixed order perturbative results as dashed lines, 
and our initial-state soft gluon resummed results as dotted lines. The asymptotic 
results are represented as dot-dashed lines. The final matched results are shown as 
solid lines. 

Fig. 6. (a) The perturbative, asymptotic, and resummed results from Fig. 5(a) are plotted 
as da/dM/dq'j_ versus q'j_. (b) The same results are plotted again but with the switching 
function Eq. (41) included as a multiplicative factor in the asymptotic and resummed 
results. 

Fig. 7. Distribution in the square of the bb pair transverse momentum. Our final matched 
results are shown by the solid line. For small qj_ we integrate the leading log resummed 
result, Eq. (|35|) , over M. The pure C(a^) result is represented as the dashed curve, the 
initial-state soft gluon resummed result as a dotted line, and the asymptotic result as 
a dot-dashed line. 

Fig. 8. Distribution in the square of the cc pair transverse momentum for proton-proton 
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collision at y/S = 38.7 GeV (corresponding to fixed target beam energy 800 GeV). 
Curves are labelled as in Fig. 7. 
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